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Abstract 

Based on the data about structure and antidiabetic activity of twenty seven vanadiunn and zinc coordination complexes 
collected from literature we developed QSAR models using the GUSAR program. These QSAR models were applied to 10 
novel vanadium coordination complexes designed in silico in order to predict their hypoglycemic action. The five most 
promising substances with predicted potent hypoglycemic action were selected for chemical synthesis and 
pharmacological evaluation. The selected coordination vanadium complexes were synthesized and tested in vitro and 
in vivo for their hypoglycemic activities and acute rat toxicity. Estimation of acute rat toxicity of these five vanadium 
complexes was performed using a freely available web-resource (http://way2drug.com/GUSAR/acutoxpredict.html). It has 
shown that the selected compounds belong to the class of moderate toxic pharmaceutical agents, according to the scale of 
Hodge and Sterner. Comparison with the predicted data has demonstrated a reasonable correspondence between the 
experimental and predicted values of hypoglycemic activity and toxicity. Bis{tert-butyl[amino(imino)methyl]carbamato}ox- 
ovanadium (IV) and sodium(2,2'-Bipyridyl)oxo-diperoxovanadate(V) octahydrate were identified as the most potent 
hypoglycemic agents among the synthesized compounds. 
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Introduction 

Vanadium is a biologically important element involved in 
numerous physiological processes. The presence of transition 
valence determines wide distribution of vanadium compounds in 
nature that are found only in chemically combined form. 
Vanadium compounds, in particular organic derivatives, are 
effective oral insulinomimetics, which inhibit lipolysis, decrease 
blood glucose levels (BGL) in animals and in clinical trials, and 
stimulate insulin secretion in experimental models of Diabetes 
MeUitus (DM) [1-5]. Although mechanisms of the insulin mimetic 
effect of vanadium complexes stiU have to be clarified, their ability 
to sensitize peripheral tissues to insulin and to reduce insulin 
resistance attracts significant attention in the context of their 
potential use for the treatment of Diabetes Type 1 (DM1), 
Diabetes Type 2 (DM2) and obesity [6,7]. 

According to the Integrity database (http: //integrity. thomson- 
pharma.com), vanadium compounds have been developed as 
drugs for treatment of Diabetes, Hypertension, Ischemia and 
Myocardial Stroke. 

The major problem limiting further clinical applications of 
inorganic vanadium compounds is a low oral absorption (1-10%). 



However, it was shown that the oral bioavailability of vanadium 
complexes with organic ligands exceeded 20^0%, which can be 
considered as acceptable level of bioavailabiKty [8] . 

Currently QSAR modeling is widely applied in rational drug 
design [9] . Despite the many examples of successful application of 
QSAR for finding and optimization of new lead compounds, only 
few QSAR studies on biological activities of vanadium complexes 
have been reported (see, e.g., [10]). Limited number of QSAR 
applications to the metal-complexes may be explained by the 
difficulties with the development of chemical descriptors, which 
reflect well the quantitative structure-activity relationships [11]. 
Recently a novel QSAR approach based on Quantitative 
Neighborhoods of Atoms (QNA) descriptors and self-consisted 
regression has been proposed; and its superiority in comparison 
with many other popular methods was shown [12]. This approach 
was realized in GUSAR program: General Unrestricted Structure- 
Activity Relationships. Since GUSAR was successfully applied for 
compounds from different chemical series and various activity/ 
property endpoints [12], it was interesting to investigate if it would 
provide the accurate prediction results for vanadium complexes 
with antidiabetic activity. Moreover, GUSAR can be used not 
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only as a tool for creation of QSAR models for the series of 
compounds with known activity values but also for estimating the 
acute rats toxicity (LD50 values) for compounds, according to 
dilTerent route of administration (oral, intravenous, intraperitoneal 
and subcutaneous) [13]. 

In this study we apply the computer program GUSAR to 
analyze the structure-activity relationships of vanadium complexes 
with antidiabetic action, predict antidiabetic activity and toxicity 
of novel designed compounds, select the most promising of them 
for synthesis and biological testing, and verify the computational 
predictions by the experiment. 

Results 

Quantitative structure-activity relationships (QSAR) 

We collected information about structures of the twenty seven 
zinc and vanadium complexes with (04), (N202), (N2S2), (02S2), 
and (S4) coordination modes and their hypoglycemic activity from 
the available sources (FUe SI). They have been reported to display 
considerable insulinomimetic elfects in in vitro experiments on 
enzyme inhibitions and blood glucose-lowering elfects. The 
hypoglycemic activity was presented as 50 percent inhibitory 
concentration (IC50). Both zinc and vanadium complexes were 
used for QSAR analysis. The collected data was randomly divided 
twenty times onto the training and test sets with a ratio of 70% and 
30%, respectively. Here we call this multiple splitting procedure a 
leave-30%-out cross-validation (L30%CV). This procedure is 
similar to n-fold external validation, but there is only one 
difference. During the n-fold validation each compound may be 
used as a test compound only once, but during the multiple 
splitting validation each compound may be used as a test 
compound several times, which depends from splitting. The 
training sets were used for creation of QSAR models and test sets 
for the assessment of external predictive accuracy. In addition to 
L30%CV we have used the procedure of selection of the test set 
compounds taking into account the activity distribution (Act.Dis.) 
using the same proportion as mentioned above: 70% for the 
training set and 30% for the test set. QSAR models were 
developed using GUSAR program (version 201 1), which based on 
Multilevel and Quantitative Neighborhoods of Atoms (MNA, 
QNA) descriptors [14,15] and self-consistent regression (SCR) 
algorithm [12,13]. 

Two hundred QNA and MNA based models were generated for 
each training set. Models that satisfied the following criteria were 
considered as acceptable for prediction: I) R^,^.st>0.5, II) Q^>0.5. 
As a result, 14 acceptable models were selected. Characteristics of 
the obtained models as well as the results of external predictivity of 
each model calculated by multiple splitting procedure (L30%CV) 
and by activity distribution procedure (L30%Act.Dis) are present- 
ed in Table 1. According to the results described in Table 1 the 
multiple splitting procedure provides more reliable assessment of 
the predictivity of models compared to the activity distribution 
procedure, which is overoptimistic. Thus, selection of models had 
been done by L30%CV. Only QNA based models were satisfied 
the described criteria. Thus, QNA descriptors are more suitable 
for modeling of complex compounds compared to MNA 
descriptors. 

Selected models were used for consensus prediction of 1 0 novel 
vanadium coordination complexes designed in silica, to estimate 
their hypoglycemic activity. Prediction results obtained using the 
consensus model for 10 novel compounds are presented in Table 2. 

Five most active compounds according to the prediction results 
described in Table 2 were selected for synthesis and hypoglycemic 
activity testing. 



Synthesis 

Based on prediction, we have selected for future synthesis five 
vanadium complexes that are expected to be more effective in 
lowering the glucose concentration in blood serum than VOSO4 
that is usually used as a reference substance [16]. Among the 
selected structures the vanadium complexes in various oxidation 
states (both neutral and anionic forms) were found, which contain 
a variety of coordination modes: oxocomplexes II, III and TV of 
vanadium(IV) - the general formula VOL2; oxocomplexes V of 
vanadium(V) with tridentate ligand, the general formula 
[V0(H20)(L)]; and oxodiperoxo complexes I of vanadium(V) - the 
general formula [V0(02) sWd" ) where n= 1. 

Complex I was formed by the interaction of sodium 
metavanadate with bidcntatc- chelating ligand 2,2'-Bipyridine 
(CioHgNg) and hydrogen peroxide: 

NaV03-l-CioH8N2-hH202 

^Na[VO(02)2(CioH8N2)] +H2O 

Formation of complex II is a two-step reaction. First, V2O5 is 
reduced to vanadyl ion [VO]^^ which is found as vanadyl sulphate 
([V(0)(H20)4]S04): 

V2O5 -I- 7H2O -h SO2 + H2SO4 
^2[V(0)(H20)4]S04 

Then, vanadyl sulphate reacts with acetyl acetone in basic 
medium: 

[V(0)(H20) J SO4 + C5H8O2 +Na2C03 

^ [V(0) (C5H702)2 + Na2S04 4H2O -H CO2 

and the overall reaction is: 

2V205-f9C5H802 

^4VO(C5H702)2 + (CH3CO)2CO + 5H2O 

Complex IV was prepared by reacting of tartaric acid with 
sodium vanadium oxide (\') under heating and further neutrali- 
zation with potassium hydroxide solution: 

V2O5 -I- 2C4H6O6 + 4KOH -I- H2O 

^K4 [(VO)2(C4H406)2] ■2H2O 60H- 

Complex III and complex V were formed according to Fig. 1 

and Fig. 2 respectively. 

The information about structure of complexes is summarized in 
Table 3. I hc chemical structures of these compounds were 
characterized by 'H-, "V-NMR spectra for vanadiumfV) 
complexes (Fig. 3A-C), IR and Raman spectra for vanadium(rV) 
and (V) complexes (Fig. 3D, E). The crystal structures of the 
vanadium complexes I, II, TV and V were confirmed by X-ray 
diffraction analysis (Fig. 4 A-E). 
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Table 2. Prediction results achieved by consensus model for 10 novel compounds. 



Compound 


Organic ligands 


IC50 Predicted 


Name of the ligand 


1 


C10N2 


0.5 


Bipyridine 


2 


C5H7O2 


0.53 


Acetyl acetone 


3 


C6H12O2N3 


0.75 


Tert-butyl[amino(imino}methyl]carbamate 


4 


C10H11NO3 


0.87 


2-(W-salicylidene)aminopropionate 




r H n 


0 93 


Tartaric acid 


VOSO4 as reference compound inorganic 
vanadium salt IC5o= 1.00±0.09 


none 


1.00 


none 


6 


C8H4CI2N2OS 


1.12 


1 ,3.4-oxadiazole-2-thiol, 5-(2,4-dichlorophenyl)- 


7 


C8H5CIN2OS 


1.13 


1 ,3,4-oxadiazole-2-thiol, 5-(4-chlorophenyl)- 


8 


C8H6N2O2S 


1.16 


3-(5-mercapto-1,3/4-oxadiazol-2-yl)- 


9 


CsHjBrNjOS 


1.16 


5-(3-bromophenyl)-1,3,4-oxadiazole-2-thlol 


10 


C8H6N2OS 


1.21 


1 ,3/4-oxadiazole-2-thiol, 5-phenyl- 



doi:l 0.1 371/journal.pone.Ol 00386.1002 



Discussion 

The developed QSAR models were applied to 10 novel 
vanadium coordination complexes designed in silico to predict 
their hypoglycemic activity. The organic molecules with donor 
atoms, preferably N and O, were used as ligands since the cations 
of vanadium (IV) and vanadium (V) form a stronger bond with 
oxygen than to nitrogen. 

It is well known that both V(IV) and V(V) complexes show 
insulin mimetic effects [6] , however it is generally considered that 
the most effective insulin-mimetic vanadium complexes are neutral 
V(IV) compounds of the general formula VO[L2]''. Several such 
VOL2 type complexes with hypoglycemic activity have already 
been patented. The bis(maltolato)oxovanadium(rV) BMOV was 
introduced into clinical tests several years ago [17]. Also, an ethyl- 
derivative of BMOV has been recently entered into clinical trials 
[18]. 

On the other hand, vanadium complexes with oxidation state + 
V electronically and structurally imitate phosphate. Recent studies 
have also demonstrated that some of pentavalent vanadium 
complexes with organic ligands are much less toxic than caffeine 
(LD50 is 192 mg/kg) and oxocompex - BMOV (LD50 is 220 mg/ 
kg) and, in addition, they are more soluble and bioavailable [17]. 

Therefore, predictions were made for various vanadium 
structures such as well known VOL2 type and for the vanadium 
(V) complexes that, in our opinion, could be also very efficient as 
insulin mimetic agents. All compounds match to the applicability 
domain of the obtained QSAR models. Five most promising 
substances with predicted potent hypoglycemic action were 
selected for chemical synthesis. 



Various approaches were applied for the synthesis of vanadium 
complexes I-V (see Methods), using as a starting material both 
inorganic vanadium(IV) and vanadium(V) compounds. 

To analyze the behavior in physiological solutions and to 
determine the structure of vanadium complexes we used various 
physico-chemical methods (IR, NMR, X-Ray). Vanadium - 51 
(nuclear spin = 7/2) is about 40% as sensitive as protons toward 
NMR observation, and therefore spectra are generally easily 
obtained (Fig. 3A). Method '''V-NMR accurately recognizes the 
specific behavior of compounds V (V) in solutions. During 
registration of '''V-NMR spectra as the external standard VOCI3 
was used. 

We used the method of '''V-NMR spectroscopy for the study of 
diamagnetic vanadium complexes I and V. The '''V-NMR 
spectrum of complex I is represented by a singlet at -662 ppm 
(Fig. 3B) typical of oxodiperoxovanadium complexes. The '''V- 
NMR spectra of complex V shows two signals in the range from — 
540 to -565 ppm (Fig. 3C) corresponding to oxovanadium 
complexes. The presence of two signals in the spectra of complex 

V can be evidence of the existence of two diastereomeric endo and 
exo forms of that complex (Fig. 5). 

Infrared spectroscopy is used mainly to prove the presence of 
fvmctional groups in the molecule. Thus, the stretching vibrations 
of the V = O in oxocomplexes of tetravalent vanadium is usually 
found at higher frequencies compared to the oxocomplexes of 
pentavalent vanadium. As an example in the IR spectra of 
complexes II the strongest bands at 980 cm ' correspond to the 

V = O stretching vibrations (Fig. 3D). In the Raman spectrum of 
complex II, the most intense lines observed at frequencies of 993 
and 466 cm"' correspond to the stretching vibrations of the V = O 
bonds and to the totally symmetric vibrations of the V-O bonds, 



II,C 



CI I, 




NH 



VOSO4 



H3C- 



-CH3 



NHt 



CII3 



Figure 1. Synthesis of oxovanadium complexes with tert-butyl [amino(imino)methyl] carbamate ligand. 

doi:10.1371/journal.pone.0100386.g001 
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HjC O 



=N 



CO 



HjC O 




J/ 



N OH 



^ VO(C5H702)2 




Figure 2. Synthesis of 2-(/V-salicylidene)aminopropionate and tKieir vanadium oxocomplex. 

doi:1 0.1 371 /journal.pone.01 00386.g002 



respectively. The medium-intensity line at 610 cm" and an 
intense line at 486 cm"' are associated with the two other among 
the three remaining of VO4 vibrations (Fig. 3E). 

The most intense bands in the IR spectrum of the synthesized 
complex I correspond to vibrations of crystal water molecules 
(3450-3200 cm" ) and to the stretching vibrations of the 
oxodiperoxovanadium moiety: the band at 926 cm"' is due to 
the V = 0 vibrations, the bands at 879 and 858 cm"' correspond 
to O-O bond vibrations. The broad band at 586 cm"' and the 
medium-intensity band at 624 cm"' are produced by the 
symmetric and antisymmetric vibrations of the VO2 fragment, 
respectively. The positions of these bands coincide with published 
data to within 4 cm"' [19]. The band at 766 cm"', which virtually 
retains its position in the spectrum of crystalline 2, 2'-Bipy, and a 
number of weak bands in the region of 1200-950 cm" can be 
attributed to ligand vibrations. The positions and intensities of 
most of the 2, 2'-Bipy bands in the IR spectrum of complex I, 
primarily in the region of 1 700-1 300 cm" , are markedly difiFerent 
from those observed for ligand-free form 2, 2'-Bipy (see 
abbreviation in Table 3), which is due to the conformational 
change that takes place in the 2, 2'-Bipy molecule upon 
complexation: the trans-arrangement of the pyridine rings in the 
free ligand is changed for the cis-arrangement in the complex. 

In the IR spectra of complex V, the most intense bands 
correspond to the vibrations of crystallization water molecules 
(3450-3200 cm"') and to the stretching vibrations of the V = O 
bonds (at 980 cm"'). In the IR spectrum of this complex, the 
ligand vibrations are primarily responsible for the band at a 
frequency of 1625 cm"', which is associated with the stretching 



vibrations of the CH = N bonds of salen ligand (see abbrevia- 
tion in Table 3). 

The crystal structures of vanadium compounds are also 
determined by using single-crystal X-ray difiFraction. The crystal 
data and refinement parameters for the structures of compounds 
are presented in the Methods. The crystal data for structures 
(CIFfiles) of complexes I, II and V have been additionally 
deposited with the Cambridge Structural Database. 

The observed values of hypoglycemic activity for vanadium 
complexes reasonably correspond to those calculated on the basis 
of QSAR models (Table 4). It was interesting that predictive value 
of IC50 (FFAs) for anionic V(V) complex I was better than the 
predictive values of IC50 (FFAs) for all neutral V(IV) structures (see 
Table 4) contrary to the statement mentioned above. 

In vitro study of the free fatty acids decrease amount released in 
adipocytes shows that FFAs level goes down in the presence of 
0.5-1 mM of vanadium complexes (see Table 4) and all of them 
were more effective than VOSO4. The most effective complexes 
according to in vitro experiment were neutral complex III which 
belongs to VOL2 type and anionic complex I which belongs to 
[VO(02)2L]" type vanadium compounds. 

Recently it was shown that there is a correlation between 
glucose level and FFAs in normal and streptozotocin-diabetic rat 
(STZ-Rats) [20,21]. It was observed that correlation coefficient of 
the liner regression (y= 177.5x+97.3) was 0.84 for a total of 64 for 
normal and STZ-Rats [20,21]. This indicates that the serum FFAs 
level along with the serum glucose level is a good index of the 
degree of diabetes mellitus. 

In order to test whether the correlation between glucose level 
(experimental data) and FFAs (predicted values) we carried out 



Table 3. Structures of synthesized vanadium complexes l-V. 





Complexes 


Formula 


Name 


Abbreviation in article 


Complex 1 {anionic} 


Na[VO(02)2(CioH8N2)]-8H20 


Sodium{2,2'-Bipyridyl) oxo-diperoxovanadate{V) octahydrate 


2,2'-Bipy 


Complex II {neutral} 


[VO(C5H702)2] 


Bis(acetylacetonato)oxo-vanadium(IV} 


acac 


Complex III {neutral} 


[V0(C6Hi202N3)2] 


Bis{tert-butyl[amino{imino} methyl] carbamato) oxovanadium (IV) 


bamc 


Complex IV (anionic) 


K4[V202{C4H406}2]-2H20 


Potassium bis{(m2-d,d-tartrato-0,0',0",0"'}-oxo-vanadium{IV)} dihydrate 


tartrate 


Complex V (neutral) 


[VO(CioHnN03}(H20)] 


(S)-[2-{N-salicylidene) amino-propionate] oxovanadium{IV} monohydrate 


salen 



doi:l 0.1 371 /journal.pone.01 00386.t003 
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in vivo study of hypoglycemic activity of our complexes I— V. In 
vivo hypoglycemic activity was determined in the model of 
adrenaline hyperglycemia caused by subcutaneous administration 
of 0.1% epinephrine hydrochloride. The results of decreasing 
blood glucose level (BGL) by vanadium compounds which were 
expressed in the percentage are summarized in Table 4. The good 
correlation (R^ = 0.91) between experimental and calculated 
values was observed (Fig. 6). 

Estimation of rodent acute toxicity is an important task in drug 
design and risk assessment of chemicals. LDjo value is one of the 
important characteristics of acute toxicity that corresponds to the 
dose causing 50% mortality within 24 hours of administration. In 
order to evaluate the toxicity of selected complexes, we used 
"GUSAR online" web service which allows predicting the acute 
rat toxicity (LD50 values) with oral administration for compounds 
based on their structures. Obtained values of LD50 for selected 
complexes and experimental data are shown in Table 5. It should 
be noted that the class of toxicity and values of LD50 for all tested 
complexes were predicted accurately. 

Conclusions 

Five coordination complexes of vanadium were synthesized and 
tested in vivo for their hypoglycemic activities and toxicity and 
were compared with the predictive data. Reasonable correspon- 
dence between the experimental and predicted values of toxicity 
and hypoglycemic activity for vanadium compounds indicates that 
GUSAR software may be successfully applied to explore the 
structure-activity relationships of metal complexes. Bis{tert- 



I 




>rv 



A 



butyl [amino(imino)methyl] carbamato } oxovanadium(IV) which 
IC50 value was almost twice less than those obtained by vanadyl 
sulfate was identified as the most potent hypoglycemic agent 
among the synthesized compounds. We have shown that it may be 
possible to develop some vanadium compounds which have 
stronger selectivity against FFAs. 

Methods 

QSAR Method 

GUSAR (version 201 1) is based on Multilevel and Quantitative 
Neighborhoods of Atoms (MNA, QNA) descriptors [13,14]. The 
calculation of QNA descriptors is based on the connectivity matrix 
(C), and also, on the standard values of ionization potential (IP) 
and electron affinity (EA) of atoms in a molecule [11,12]. 

For any given atom i, the QNA descriptors are calculated as 
follows: 



Qi = B,Y{Exp(-\c)) BUk, 



with Ak = -(IPk+EAk), Bt = (IPk-EAk)-K 




B C 




I I I I I I I I MX) :iw 2IXtJ iMI luuu :<ju 

40U CiUU SUU lUUU 12UU l-MJO 16UU liiUU ItJUO v.cm-' 

V. cur' 



D E 

Figure 3. SIV-NWIR spectra for vanadium (V). The wide spectral dispertion of the signals as characteristic of vanadium NMR spectra (A) complex 
I (B); complex V (C). IR (D) and Raman (E) spectra for vanadium complex II. 
doi:1 0.1 371/journal.pone.01 00386.g003 
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E 

Figure 4. Molecular structure and the atomic numbering for complex 1(C), complex 11(A), complex IV(D, E) and complex V(B). The 

ellipsoids of thermal vibrations are shown at the 50% probability level. 
doi:1 0.1 371/journal.pone.01 00386.g004 



The P and values are calculated for aU atoms of the molecule. 
Two-dimensional Chebyshev polynomials are used for approxi- 
mating the functions P and Q. Thus, the independent regression 
variables are calculated as average values of particular two- 
dimensional Chebyshev polynomials of P and values for the 
atoms in a molecule. 

In addition, GUSAR allows the creation of QSAR models 
based on predicted biological activity profiles of compounds. This 
is done by using the PASS algorithm on each compound's 



representation as a list of MNA descriptors for the prediction of 
the compound's biological activity profile [22,23]. The PASS 
program version of 10.1 predicts 4130 types of biological activity 
with a mean prediction accuracy of about 95%. The list of 
predicted biological activities includes pharmacotherapeutic ef- 
fects, mechanisms of action, adverse and toxic effects, metabolic 
terms, susceptibility to transporter proteins and activities related to 
gene expression. The results of the PASS procedure are output as 
a list of the difference between the probability, for each biological 
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endo-\ exo-V 



Figure 5. Structures of two diastereomeric forms of complex V, 
R = CH3, R' = 0. 

doi:1 0.1 371 /joumal.pone.01 00386.g005 

activity, of the compound to be active (PJ or to be inactive (P;). For 
obtaining the different QSAR models in GUSAR, subsets of these 
Pa-P; values were randomly selected from the total list of predicted 
biological activities as input independent variables for the 
regression analysis. The size of these subsets depends on the 
number of compounds in the training set. If the number of 
compounds in the training set falls between 10 and 30, then the 
number of initial variables is closer to the number of compounds in 
the training set. 

GUSAR (version 2011) also allows the calculation of whole- 
molecule descriptors: topological length, topological volume and 
lipophUicity. These descriptors are used in combination with the 
QNA and MNA descriptors. 

For generation of the QSAR models, GUSAR uses a self- 
consistent regression (SCR) algorithm. SCR is based on the 
regularized least-squares method. Unlike stepwise regression and 
other methods of combinatorial search, the initial SCR model 
includes all regressors. The basic purpose of the SCR method is to 
remove the variables that poorly describe the modeled value but to 
retain the set of variables correcdy representing the existing 
relationship. The number of final variables in the QSAR equation 
selected after the SCR procedure is typically significandy lower 
than the number of the initial variables. The details of the 
algorithms for descriptor calculation and the self-consistent 
regression methods have been described previously [22,23]. 

For assessment of the applicability domain of the obtained 
models GUSAR uses the three different approaches: similarity, 
leverage, and accuracy assessment. The detailed description of 
these approaches is presented here [24]. 

The QSAR models obtained using both MNA and QNA 
descriptors are integrated into the consensus model which is 
further used for prediction. The final predicted value is estimated 
by taking into account a weighted average of the predicted values 
from each obtained QSAR model (QSAR models provide the 
predictions that are within the respective applicability domains). 
The predicted value obtained from each developed model is 
weighted on similarity value calculated during the evaluation of 
applicability domain. 

Chemistry Method 

General methods for preparation of the vanadium 
derivatives. AH reagents and solvents were used without further 
purification or drying. All the reagents (including reference 
material VOSO4) were purchased from Vekton (Russia). All 
reactions were proceeded in thermostat reactor under argon 
atmosphere. Beforehand, the solvents were dried and distilled 
under nitrogen using standard methods. Elementary analyses of H, 
C and N were performed with a Carlo Erba microanalyzer. 'H- 
NMR spectra at 400MHz were recorded on Varian XR-400 
instrument, '^'V-NMR spectra at 52.6MHz were recorded on 
Tesia Bruker AL-200 instrument, and mass spectra were obtained 



on a Finnigan MAT 1 12S spectrometer. IR spectra were obtained 
on a Specord M80 (Carl Zeiss) spectrophotometer within the 
spectral range of 400-4000 cm~' under resolution of 2.5-4 cm~'. 
Spectra of the solutions in CH2CI2 or CHCI3 were shown to 
represent pure solvents; the films were obtained on a KBr 
substrate by evaporation of the solvent. X-Ray diffraction study of 
the single crystal was performed on an Enraf-Nonius CAD4 
automated diffractometer. The unit cell parameters were deter- 
mined and refined using 25 reflections in the range of 15°^6^ 
16°(XMoKa radiation graphite monochromator). 

Synthesis of compound I - Na[VO(02)2(bipy)] -SHzO. Com- 
plex I was prepared following the procedure reported in [16]. A 
solution of 2,2'-pipydidine (0.8 g) in ethanol (10 ml) was added to a 
solution of sodium metavanadate in cooled 20% H2O2 (20 ml). 
After 2 to 3 min, ethanol (20 ml) was added with stirring and the 
solution was allowed to crystallize at 5°C. The yellow crystals thus 
formed were filtered off, washed once with ethanol, and dried. 
According to the data of elemental analysis, the resulting compound 
I had the composition Na[VO(02)2(bipy)]-8H20. The crystals are 
triclinic: ffl = 7.200(4) A, 6= 11.245(4) A, c= 13.703(5) A, 
a = 111.66(2)°, (5 = 90.86(2)°, y = 88.01(2)°, Z = 2, F= 1030.5(8) A, 
space group 1, Pcakd = 1 -464 g/cm'', |j,mo =0.564 mm '. 'H 
NMR, "^'V NMR spectra, IR spectra and X-ray data of the 
synthesized compound I have been described in details early by 
Fedorova at al [25]. 

Synthesis of compound II - [VO(C5H702)2]. The oxova- 
nadium complex \yO[acac)-2\ (II) was synthesized according to a 
modified procedure described in [1 7] . A mixture of V2O5 (20 g, 
0. 1 1 mol), distilled water (50 ml), concentrated sulfuric acid 
(30 ml), and ethanol (100 ml) was boiled with stirring for 
30 min. The dark blue solution obtained was filtered, and freshly 
distilled acetylacetone (50 ml, 0.49 mol) was added to the filtrate. 
Then, a solution of sodium carbonate (80 g of Na2C03 in 500 ml 
of water) in water (50 ml) was slowly added to the reaction 
mixture. The blue product was filtered off, washed with water, and 
dried in air. The yield was 50 g. Blue crystals were obtained by 
recrystaUization from chloroform. The crystals were insoluble in 
water and soluble in CH2CI2 and dimethyl sulfoxide. The crystal 
structure of compound II was determined using single-crystal X- 
ray diffraction. Crystals of compound II are triclinic, a = 
7.4997(19) A, 6 = 8.2015(15) A, c= 11.339(3) A, a = 91.37(2) °, 
P= 110.36(2)°, y= 113.33(2) °, Z = 2, and space group PI. 
NMR, ^'V NMR spectra, IR spectra and X-ray data of the 
synthesized compound II have been described in details early by 
Fedorova at al [26]. 

Synthesis of compound III - [VO(C6Hi202N3)2]. Complex 
III was prepared following the procedure reported in [18]. A 
solution of vanadyl sulfate VOSO4.2H2O in water was slowly 
added to the solution of tert-butyl carbamimidoylcarbamate in 
molar ratio 1:2,2 with continuous stirring. Thereafter, the solution 
acquired a dark green color. In a matter of minutes, a heavy 
precipitate was formed. After stirring for 30 min, the substance 
was filtered off and washed with distilled water, 50% ethanol 
solution, and ethyl ether. The compound synthesized was dried in 
a vacuum desiccator at a temperature of 40°C. The yield of the 
product was approximately equal to 75% with respect to 
VOS04-2H20. The final product was recrystaUized from 
methanol. 

Synthesis of compound IV - K4[V202{C!4H406)2] • 
2H2O. Complex rV was prepared following the procedure 
reported in [19]. A mixture of V2O5 distilled water and water 
solution of Z)-tartaric acid was boiled with stirring for 30 min. 
Then, a solution of potassium hydroxide in water was slowly 
added to the reaction mixture. The yield was 45%. The crystal 
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Table 4. Predictive values of IC50 
decreasing blood glucose level. 


in mlVl (FFAs) for vanadium compounds and 


in vivo observations IC50 in 


mM (FFAs) and 




Compounds 


V0S04 


Complex 1 


Complex II 


Complex III 


Complex IV 


Complex V 


1^50 predicted 




0.50 


0.75 


0.53 


0.93 


0.87 


ICso e^p (±S.D.) 


1.0 ±0.09 


0.72±0.09 


0.95 ±0.09 


0.54 ±0.07 


0.95±0.10 


0.81 ±0.06 


decreasing blood glucose level 


16% 


50% 


33% 


42% 


11% 


26% 



doi:1 0.1 371 /journal.pone.Ol 00386.t004 



Structures of compound IV was determined using single-crystal X- 
ray difTraction. Complex IV with square-pyramidal coordination 
geometry, space group P4(3)22, a = 7.9345 A, b = 7.9345 A, 
c = 30.5035 A, F= 1920.39(22) A^ Z = 4, R„ = 3.11%. 'H 
NMR, "''V NMR spectra, IR spectra and X-ray data of the 
synthesized compound IV have been described in details earlier by 
Fedorova at al [27]. 

Synthesis of compound V - [VO(CioH703)(H20)]. Oxo- 
vanadium complex V was synthesized according to a modified 
procedure described earlier in [17]. The ligand in the complex is a 
bivalent tridentate SchifPs base, which was produced by the 
reaction of salicylaldehyde {Sal) with ot-amino acid. L-alanine 
amino acid (0. 1 mol) and sodium acetate (0.2 mol) were dissolved 
in distilled water (200 ml). In order to dissolve the amino acid 
completely, the solution was heated and filtered. A solution of 
salicylaldehyde (0.1 mol) in ethanol (250 ml) was added to the 
filtrate. A solution of vanadyl sulfate VOS04-2H20 (0.085 mol) in 
water (80 ml) was slowly added to the above solution with 
continuous stirring. Thereafter, the solution acquired a dark 
brown color. In a matter of minutes, a heavy precipitate was 
formed. After stirring for 30 min, the substance was filtered off and 
washed with distilled water, 50% ethanol solution, and ethyl ether. 
The compound synthesized was dried in a vacuum desiccator at a 
temperature of 40°C. The yield of the product was approximately 
equal to 80% with respect to VOS04-2H20. The final product 
was recrystallized from metha similarity nol. The blue crystals thus 
prepared were insoluble in water, acetone, ether, and benzene and 
soluble in methanol, pyridine, methylene chloride, and chloro- 
form. The melting temperature of the crystals was estimated as T,-,-, 
~ 250°C. Crystals of compound V are monoclinic, a = 8.5106(16) 
A, & = 7.373(2) A, f = 9.1941(16) A, (5=101.88(1), Z = 2, and 
space group P2,. 'H NMR, "''V NMR spectra, IR spectra and X- 
ray data of the synthesized compound V have been described in 
details early by Fedorova at al [26]. 



Pharmacology Method 

All procedures were performed according to institutional 
guidelines for animal experimentation; protocol was submitted 
and approved by the Institutional Committee for Laboratory 
Animal Use and Care of the Saint-Petersburg Chemical Pharma- 
ceutical Academy. Adult male Wistar rats and mongrel mice 
obtained from the ^Rappolovo' lab animals breeding farm near St. 
Petersburg, belonged to Russian Academy of Medical Science 
were used. All experiments were approved by the Animal Ethics 
Committee of the Saint-Petersburg Chemical Pharmaceutical 
Academy. 

Insulin-mimetic action of vanadium complexes. All rats 
(250-300 g) and mice (22-25 g) were housed in a temperature 
controlled room (24°C± 1°C) and they had received standard food 
pellets and tap water until they weighed. At the end of the study 
(24 h after CIN induction), all animal test subjects were 
euthanized by inhalation of ether anesthesia. Rat adipocytes were 
prepared from the fat pads of male Wistar rats by coUagenase 
digestion according to the method of Rodbell [28]. Freshly 
prepared adipocytes were incubated (37 'C, 30 minj with solutions 
of the vanadium compounds at three concentrations (0.1, 0.5 and 
1 mM in isotonic saline and 5 mM glucose), followed by 3 hours 
of incubation with epinephrine. The inhibitory effect was 
determined in terms of the decrease of the amount free fatty 
acids released in adipocytes. All assays were performed in 
duplicate or triphcate. Table 4 summarizes the IC50 values, 
which are defined as the concentrations c(V) at which 50% 
inhibition of FFA release takes place. Additionally, five coordina- 
tion complexes of vanadium and reference material VOSO4 have 
been tested in vivo for their hypoglycemic activities. Investigations 
were carried out on mongrel white mice. 

Hyperglycemic activity. Determination of hypoglycemic 
activity of 5 complexes was performed on 75 mice. Adrenaline 
hyperglycemia was caused by subcutaneous administration of 
0.1% epinephrine hydrochloride (AGH) at a dose of 1 mg/kg. 
With a single run all the compounds were administered orally 1 
hour prior to the AGH, at 1 / 1 0 of the LD50. In each experimental 
group there were 3 animals and 3 replicates were performed. The 



Table 5. Predicted and observed values LDso 


of vanadium compounds. 






# Compounds 


LDsopred/Class in AD 


LDsoexp/Class In AD 


1 VOSO4 


518.6/3 


448(356-540)/3 


2 Complex 1 


279.0/3 


1 08(95-1 22)/3 


3 Complex II 


251.0/3 


250(21 0-280)/3 


4 Complex III 


211.0/3 


190(170-210)/3 


5 Complex IV 


357.1/3 


261(210-312)/3 


6 Complex V 


265.1/3 


245(21 0-280)/3 


doi:1 0.1 371 /journal.pone.Ol 00386.t005 
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% decreasing blood glucose level 



Figure 6. Correlation between experimental and calculated 
values for vanadium compounds. 

doi:l 0.1 371/journal.pone.01 00386.g006 

content of glucose in whole blood was determined at 0.5, 1 and 2 
hours after administration of AGH or glucose using portable 
glucose monitoring systems, OneTouch Ultra test strip and 
OneTouch Ultra and Clever check TD4209 and strips Clever 
check TD4209. Oral administration of all compounds elicited a 
progressive reduction in plasma glucose over 2 hours in mice. The 
results (as % decreasing blood glucose level) are summarized in 
Table 4. 

Acute toxicity studies. Acute toxicity studies were carried 
out by rapid method of Prozorovskii [29] . The method consists of 
se(]ucntial administration of multiple doses with logarithmic step, 
resulting in a set number of responses and in determining the value 
of LD50, its range of variation of the average error using 4 groups 
of 3 observations in each group. The following criteria were used 
for the evaluation of acute toxicity: the number of dead animals, 
the timing of their death, LD-jq through the night, the picture of 
intoxication, behavioral change and autopsy results of the survived 



animals after 14 days. Chnical examination was performed for 
each animal during the first hour after administration of the 
substance, and weekly thereafter. Clinical examination included a 
detailed inspection of the animal in cage, in hands and in the open 
area. The manifestation and severity of intoxication were noted. 
Examination of animals in order to identify mortality or signs of 
variations in health status was conducted on a daily basis. 
Integrated indicators were assessed weekly during the entire period 
of observation. Death of animals was recorded daily during entire 
study period. General condition and behavior of animals as well as 
type of toxic manifestations, including their de\-cloj)m(;nt and 
termination were recorded daily during the procedures for 
monitoring and inspection of individual animals. 

Statistical analysis. For all recorded indicators the following 
parameters of descriptive statistics were introduced: the arithmetic 
mean, the standard deviation and the standard error of the mean. 
The following statistical tools were used as methods for compar- 
ison: U- Mann-Whitney and Student's tests; the Friedman and 
Wilcoxon test for repeated measurements; the U- Mann-Whitney 
for subsequent pairwise (post-hoc) comparisons. The level of 
comparison significance in this study was set up at 5%. Results 
were considered statistically significant if the p-value for the test 
was less than or equal to 0.05. 
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